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The SciBooNE Collaboration reports a measurement of inclusive charged current interactions of 
muon neutrinos on carbon with an average energy of 0.8 GeV using the Fermilab Booster Neu- 
trino Beam. We compare our measurement with two neutrino interaction simulations: NEUT and 
NUANCE. The charged current interaction rates (product of flux and cross section) are extracted 
by fitting the muon kinematics, with a precision of 6-15% for the energy dependent and 3% for 
the energy integrated analyses. We also extract CC inclusive interaction cross sections from the 
observed rates, with a precision of 10-30% for the energy dependent and 8% for the energy inte- 
grated analyses. This is the first measurement of the CC inclusive cross section on carbon around 
1 GeV. These results can be used to convert previous SciBooNE cross section ratio measurements 
to absolute cross section values. 

PACS numbers: 13.15.+g, 25.30.Pt 
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The neutrino charged current (CC) interaction is 
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nucleus scattering in the few-GeV region. It is impor- 
tant for neutrino oscillation measurements because the 
neutral current (NC) interaction is flavor-blind. How- 
ever, this interaction mode is poorly understood because 
of large neutrino flux and cross section uncertainties. 

In the energy region above ~3 GeV, NOMAD [l| and 
MINOS recently reported precise CC inclusive in- 
teraction cross section measurements. However, in the 
~1 GeV region, all inclusive CC measurements have been 
made on deuterium targets using bubble chambers [H, 0] ■ 
In this energy region, the nuclear effects of the neutrino 
target material (from Fermi motion and the nuclear po- 
tential) are significant. Therefore, the cross section on 
deuterium targets is not directly applicable to the heavier 
nuclear target materials used in the recent accelerator- 
based neutrino experiments in this energy region, e.g. 
SciBooNE, MiniBooNE [| and T2K @|. Furthermore, 
the results in Ref. [3] use CC quasi elastic (QE) inter- 
actions to normalize the absolute neutrino flux, which 
introduces additional ambiguity from the choice of CC- 
QE interaction model parameters. 

Because of the poor knowledge of the cross section, 
there are multiple neutrino interaction simulators used 
to predict different cross sections and kinematics of fi- 
nal state particles. Among these, NEUT @, Hi and NU- 
ANCE [9( are commonly used for the recent neutrino os- 
cillation and interaction measurements. NEUT is used in 
the Kamiokande 0, Super-Kamiokande [ll|, K2K Q3, 
and T2K experiments, while NUANCE is used in Mini- 
BooNE and as a check of simulations by multiple experi- 
ments. They are both tuned to describe the experimental 
data, however, they have not yet been precisely compared 
with each other in a single experiment. 

In this paper, we report on the comparison of CC 
candidate events observed in SciBooNE with predictions 
based on the NEUT and NUANCE neutrino interaction 
simulators, as well as on the measured CC interaction 
rate (product of flux and cross section) as a function of 
neutrino energy extracted from these comparisons. We 
also report on the first measurement of the CC inclusive 
interaction cross section on carbon in the 1 GeV region, 
which is relevant to ongoing and future neutrino oscilla- 
tion experiments. 

In addition to this general purpose, this analysis is 
also motivated by two direct applications to our mea- 
surements. The first is to provide a constraint on the 
product of flux and cross-section for the forthcoming 
SciBooNE and MiniBooNE joint disappearance anal- 
ysis The MiniBooNE collaboration recently per- 
formed a search for and disappearance using only 
MiniBooNE data, in which the largest uncertainty stems 
from the flux and cross section uncertainties [14j. As 
shown in Fig. [TJ the SciBooNE detector is located 440 m 
upstream of the MiniBooNE detector, sharing the same 
neutrino beam. In addition, the neutrino target mate- 
rials are both essentially carbon; polystyrene (CsHg) for 
SciBooNE, and mineral oil (CH 2 ) for MiniBooNE. There- 
fore, most of the flux and cross section uncertainties can- 



cel when we compare the data from the two experiments. 
Then, a disappearance search with higher sensitivity 
becomes possible [13| . 
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FIG. 1. (color online). Schematic overview of the Booster 
Neutrino Beamline and the location of the SciBooNE and the 
MiniBooNE detectors, (not to scale). 

The second motivation is to provide an absolute nor- 
malization factor for SciBooNE's previous cross section 
ratio measurements. The SciBooNE collaboration re- 
cently measured relative cross sections of charged current 
coherent pion production [l5| and neutral current neutral 
pion production [Tril [l7j , which were both normalized to 
CC inclusive interactions. This paper provides an abso- 
lute CC interaction cross section with the same definition 
used in the previous analyses, so that the measured ratios 
can be converted to absolute cross sections. 

The remainder of this paper is organized as follows: 
Section llll summarizes the SciBooNE experiment, includ- 
ing the Booster Neutrino Beamline (BNB) and the Sci- 
BooNE detector. The simulation of neutrino interactions 
with nuclei is described in Section IIIII Section IIVI de- 
scribes the methods used to reconstruct CC interactions 
and defines the sub-samples used. The analysis of the en- 
ergy dependent neutrino interaction rate is described in 
Section [V] Finally, the results of the CC interaction rate 
analysis and the extracted CC absolute cross section are 
presented in Section IVII The final conclusions are given 
in Section Ivni 



II. SCIBOONE EXPERIMENT 

A. Neutrino Beam 

SciBooNE detected neutrinos produced by the Fermi- 
lab BNB. The same BNB beam is also serving the Mini- 
BooNE experiment. The BNB uses protons accelerated 
to 8 GeV kinetic energy by the Fermilab Booster syn- 
chrotron. Beam properties are monitored on a spill-by- 
spill basis, and at various locations along the BNB line. 
Transverse and directional alignment of the beam, beam 
width and angular divergence, beam intensity and losses 
along the BNB, arc measured and used in the data qual- 
ity selection [TEj. Protons strike a 71.1 cm long beryllium 
target, producing a secondary beam of hadrons, mainly 
pions with a small fraction of kaons. A cylindrical horn 
electromagnet made of aluminum surrounds the beryl- 
lium target to sign-select and focus the secondary beam. 
For the data set used in this measurement, the horn po- 
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larity was set to neutrino mode, focusing particles with 
positive electric charge. The neutrino beam is mostly 
produced in the 50 m long decay region, mainly from 
7r + — > in-flight decays. 

The analysis presented here uses the full neutrino data 
set: 0.99xf0 20 protons on target (POT) collected be- 
tween October 2007 and April 2008. 



B. Neutrino Flux Prediction 

Predictions for the BNB neutrino flux illuminating the 
SciBooNE detector are obtained via a GEANT4 [H] sim- 
ulation of the beamline. Hadronic interactions in the 
beryllium target are carefully modeled. For ir + produc- 
tion, a parametrization based on HARP [l9| and BNL 
E910 [20j data are used. The total, inelastic and quasi- 
elastic hadronic cross sections of protons, neutrons, 7r + s 
and 7T~s with Be or Al are treated via custom cross sec- 
tion models [2l| . Other hadronic and all electromagnetic 
processes of importance to neutrino production are de- 
scribed by standard GEANT4 models. For a detailed de- 
scription of the BNB simulation code, see Ref. [2l[ . A to- 
tal neutrino flux per proton on target of 2.2 x 10~ 8 cm~ 2 
is expected at the SciBooNE detector location for the 
neutrino running mode. The flux is dominated by muon 
neutrinos (93%). 
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FIG. 2. (color online). (Top) flux prediction at the Sci- 
BooNE detector as a function of neutrino energy E v . The to- 
tal flux and contributions from 7r + and K + decays are shown. 
(Bottom) Fractional uncertainty of the flux prediction due 
to 7r + and K + production from the p-Be interaction. Ad- 
ditional uncertainties from POT, hadron interactions in the 
target, and the horn magnetic field are not shown. 

Figure [5] shows the flux from 7r + and K + de- 



cays, and their fractional uncertainties. The energy 
spectrum peaks at ~ 0.6 GeV, and extends up to 2- 
3 GeV. The mean energy of the flux is expected to 
be 0.76 GeV. The largest error on the predicted neutrino 
flux results from the uncertainty of pion production in the 
initial p-Be process in the target. The simulation predicts 
that 96.7% of muon neutrinos in the BNB are produced 
via 7r + decay. The uncertainty in 7r + production is deter- 
mined from spline fits to the HARP 7r + double differential 
cross section data [2l[. The HARP data used are those 
from a thin (5% interaction length) beryllium target run 
19]. While the HARP data provide a valuable constraint 
on the BNB flux prediction, additional uncertainties re- 
sulting from thick target effects (secondary re-scattering 
of protons and pions) are included through the BNB flux 
simulation. The resulting 7r + production uncertainty is 
5% at the peak of the flux distribution and increases 
significantly at high and low neutrino energies. 

The flux from K + decay is dominant for E u > 2.3 GeV. 
Since no published data exist for K + production at the 
BNB primary proton beam energy, we employ the Feyn- 
man scaling hypothesis to relate K + production measure- 
ments at different proton beam energies to the expected 
production at the BNB proton beam energy [2l| . The 
errors of the Feynman scaling parameters obtained from 
these measurements are accounted for in the systematic 
error. Other major contributions to the flux error in- 
clude uncertainties on hadron interactions in the target 
and simulation of the the horn magnetic field, which both 
contribute to shape and normalization uncertainties, as 
well as the measurement of the number of POT, which 
is a pure normalization uncertainty. All flux errors are 
modeled through variations in the simulation and result 
in a total error of « 7% at the peak of the flux. 



C. SciBooNE Detector 

The SciBooNE detector is located 100 m downstream 
from the beryllium target on the beam axis, as shown 
in Fig. [TJ The detector comprises three sub-detectors: 
a fully active and finely segmented scintillator tracker 
(SciBar), an electromagnetic calorimeter (EC), and a 
muon range detector (MRD). SciBar is the primary neu- 
trino target for this analysis. In this analysis, we use the 
data from SciBar and the MRD, and the data from the 
EC are not used. 

SciBooNE uses a right-handed Cartesian coordinate 
system in which the z axis is the beam direction and 
the y axis is the vertical upward direction. The origin 
is located on the most upstream surface of SciBar in the 
z dimension, and at the center of the SciBar scintillator 
plane in the x and y dimensions. Since each sub-detector 
is read out both vertically and horizontally, two views 
are defined: top (x vs. z projection) and side (y vs. z 
projection). 

The SciBar detector [22j was positioned upstream of 
the other sub-detectors. It consists of 14,336 extruded 
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plastic scintillator strips. Each strip has a dimension of 
1.3 x 2.5 x 300 cm 3 . The scintillators are arranged ver- 
tically and horizontally to construct a 3 x 3 x 1.7 m 3 
volume with a total mass of 15 tons. The dominant com- 
ponent of the SciBar detector is polystyrene (CsHg). We 
measured the density of the scintillator by sampling 10% 
of the strips before the installation [23| . The uncertainty 
of the total detector mass is estimated to be 1%, includ- 
ing the effect of epoxy resin used to glue the strips. 

Each strip is read out by a wavelength shifting (WLS) 
fiber attached to a 64-channel multi-anode photomulti- 
plier (MA-PMT) . Charge information is recorded for each 
channel, while timing information is recorded in groups 
of 32 channels by taking the logical OR with multi-hit 
TDC modules [24j . The timing resolution for minimum- 
ionizing particles, evaluated with cosmic ray data, is 
1.6 ns. The average light yield for minimum-ionizing par- 
ticles is approximately 20 photo-electrons per 1.3 cm path 
length, and the typical pedestal width is below 0.3 photo- 
electron. The hit finding efficiency, evaluated with cosmic 
ray data, is more than 99.8%. The minimum length of a 
reconstructable track is approximately 8 cm (three layers 
hit in each view). The track finding efficiency for single 
tracks of 10 cm or longer is more than 99%. 

The MRD is located downstream of SciBar and the EC, 
and is designed to measure the momentum of muons pro- 
duced by CC neutrino interactions. It comprises 12 iron 
plates with thickness 5 cm sandwiched between planes 
of 6 mm thick scintillation counters; there are 13 alter- 
nating horizontal and vertical planes read out via 362 
individual 2 inch PMTs. Each iron plate measured 274 
x 305 cm 2 . The MRD measures the momentum of muons 
up to 1.2 GeV/c using the observed muon range. Charge 
and timing information from each PMT are recorded. 
The average hit finding efficiency is 99%. 



D. Detector Simulation 

A GEANT4 framework is used for the detector simula- 
tion. The detector simulation includes a detailed geomet- 
ric model of the detector, including the detector frame 
and experimental hall and soil. 

The Bertini cascade model within GEANT4 [H| is 
used to simulate the interactions of hadronic particles 
with detector materials. A 10% difference of the total 
pion-carbon cross section is seen for higher energy pi- 
ons between the GEANT4 simulation and external mea- 
surements. Hence, we set ±10% systematic uncertainty 
for the pion-carbon cross section. The systematic un- 
certainty of the energy deposit per unit length is esti- 
mated to be 3% in SciBar and MRD, and 10% in EC, 
which is dominated by the differences among various 
calculations of the range to energy conversion. Birk's 
constant for the SciBar scintillator was measured to be 
0.0208 ± 0.0023 cm/MeV 

The crosstalk of the MA-PMT is measured to be 
3.15±0.4% for adjacent channels. The single photo- 



electron resolution of the MA-PMT is set to 50% in the 
simulation, to reproduce the observed dE/dx distribu- 
tion of cosmic muons. The absolute error is estimated 
to be ±20%. A more detailed description of the detector 
simulation is given in (l5j . 



III. NEUTRINO INTERACTION SIMULATION 
A. Overview 

We simulate neutrino interactions with carbon and hy- 
drogen in the SciBar detector with the NEUT and NU- 
ANCE program libraries. We produced separate Monte 
Carlo samples with NEUT and NUANCE, and compared 
the two to the SciBooNE neutrino data. 

The nuclear targets handled in NEUT are protons, car- 
bon, oxygen, and iron. The energy of neutrinos handled 
by the simulation ranges from 100 MeV to 100 TeV. 
The types of neutrino interactions simulated in both 
NC and CC are : elastic and quasi-elastic scattering 
(uN — > IN'), single meson production (vN — > IN'm), 
single gamma production (vN — > iN'j), coherent 7r pro- 
duction (> 12 C(or 56 Fe) ->• Zir 12 C(or 56 Fe)), and deep in- 
elastic scattering (vN — > (N'hadrons), where TV and 
N' are the nucleons (proton or neutron), i is the lep- 
ton (electron, muon or neutrino), and m is the meson. 
In nuclei, interactions of the mesons and hadrons with 
the nuclear medium are simulated following the neutrino 
interactions. 

The types and models of neutrino interactions in NU- 
ANCE are similar to those in NEUT but with different 
cross section parameter settings in some cases and a com- 
pletely independent treatment of meson and hadron rc- 
intcractions in the nuclear medium. 

In addition to the neutrino interactions inside SciBar, 
we also simulate interactions in the EC/MRD and the 
surrounding materials (the walls of the detector hall and 
soil) using NEUT. We assign a 20% normalization un- 
certainty to the interaction cross sections for both the 
EC/MRD and surrounding material relative to the pre- 
dictions for SciBar. 



B. Neutrino Interaction Model 

1. Quasi-elastic scattering 

The dominant interaction in the SciBooNE neutrino 
energy range is CC-QE scattering, which is implemented 
using the Smith and Moniz model [13] ■ The nucleons are 
treated as quasi-free particles and the Fermi motion of 
nucleons along with the Pauli exclusion principle is taken 
into account. The Fermi surface momentum (j>f) for car- 
bon is set to 217(220) MeV/c and the nuclear potential 
(E B ) is set to 25(34) MeV/c in NEUT (NUANCE) , as 
extracted from electron scattering data [28j . The default 
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binding energy in NUANCE is somewhat higher because 
it additionally accounts for neutrino vs. electron scatter- 
ing differences [13 ■ The systematic errors for pp and Eb 
are set to ± 30 MeV/c and ± 9 McV/c, respectively, for 
both NEUT and NUANCE. 

For the vector form factor, NEUT uses a dipole form 
with a vector mass of 0.84 GeV/c 2 , while NUANCE uses 
the BBA-2003 form factor [30]. A dipole form is used 
for the axial form factor with an adjustable axial mass, 
M% E , for both NEUT and NUANCE. In NUANCE, an 
empirical Pauli-blocking parameter, k, is introduced (3lj | 
to better describe the MiniBooNE quasi-elastic data at 
low momentum transfer. When k > 1, the phase space 
of nucleons susceptible to Pauli-blocking is increased and 
hence the cross section at low momentum transfer is sup- 
pressed. 

The values of M% E = 1.21 GeV/c 2 and k = 1.000 
(i.e. no additional Pauli blocking adjustment) are used 
in NEUT, and M% E = 1.23 GeV/c 2 and k = 1.022 
are used in NUANCE |29j. A systematic uncertainty of 
± 0.22 GeV is assigned to M® E to span the difference 
between the value used and the global fit from previous 
The difference between k 



measurements |32j|. The difference between k = 1.000 
and k = 1.022 is also assigned as systematic uncertainty. 

The same Fermi momentum distribution, nuclear po- 
tential are used in all other neutrino-nucleus interactions 
except for coherent 7r production. 



2. Meson production via baryon resonances 

The second most frequent interaction in SciBooNE is 
the resonant production of single pion, kaon, and eta 
mesons as described by the model of Rein and Sehgal 

(RS) m. 

The RS model assumes an intermediate baryon reso- 
nance, N*, in the reaction of vN — > £N*,N* — > N'm. 
All intermediate baryon resonances with mass less than 
2 GeV/c 2 are included. Baryon resonances with mass 
greater than 2 GeV/c 2 are simulated as deep inelastic 
scattering. A re-interactions (AN — > NN) which do not 
lead to a mesonic final state are also simulated. This 
re-interaction probability is assumed to be 0.2 ± 0.2 for 
all A resonances in NEUT and 0.1 ± 0.1 (0.2 ± 0.2) for 
resonances in NUANCE. 

To determine the angular distribution of final state pi- 
ons, the RS method (34| is used for the P33(1232) res- 
onance in both NEUT and NUANCE. For other reso- 
nances, the directional distribution of the generated pion 
is chosen to be isotropic in the resonance rest frame. 

The axial-vector form factors are formalized to be 



dipole with M 1 / = 1.21 GeV/c 2 for NEUT and Aff = 
1.10 GeV/c 2 for NUANCE, with an uncertainty of 
0.28 GeV/c 2 in both cases. 

An additional uncertainty is assigned to account for 
the observed Q 2 disagreement between the SciBooNE CC 
l7r-enriched data samples and NEUT [ll| . A similar dis- 
agreement is also observed for NUANCE. The size of the 



uncertainty is determined by re-weighting CC resonant 
pion events as a function of true Q 2 such that they match 
the observed distribution in the SciBooNE data, and us- 
ing the difference between the re-weighted distribution 
and the central value as the uncertainty. 

Resonance decays leading to multi-pion final states are 
also included in the model and are simulated assuming 
= 1-30 ± 0.52 GeV/c 2 in NUANCE. This value 
of M j 4 > is chosen strictly to ensure that the total CC 
cross section prediction reproduces previous experimen- 
tal data. In NEUT, multi-pion production is simulated as 
deep inelastic scattering as described later in Sec. lIIIB~4l 
and the RS model is not used. The size of the system- 
atic uncertainty is estimated based on variations 
for NUANCE, and the same size error is assumed also 
for the NEUT prediction. 



3. Coherent pion production 

Coherent pion production is a neutrino interaction 
with a nucleus which remains intact, releasing one pion 
with the same charge as the incoming weak current. Be- 
cause of the small momentum transfer to the target nu- 
cleus, the outgoing pion tends to be emitted in the for- 
ward direction, closely following the incoming neutrino 
direction. The formalism developed by Rein and Seh- 
gal [H, 36] is used to simulate such interactions. The ax- 
ial vector mass, M%> h , is set to 1.0±0.28 GeV/c 2 in both 
NEUT and NUANCE. In NEUT, the total and inelastic 
pion-nucleon cross sections from the original Rein-Sehgal 
publication are employed (35l.l36T|. In NUANCE, they are 
obtained from fits to PDG data [37] and implemented as 
a function of pion energy. Additionally, the NC and CC 
coherent pion production cross section predictions in NU- 
ANCE are rescaled by a factor of 0.65 to better match 
the measured rate of NC coherent ir° production as mea- 
sured in MiniBooNE f38j . 



4- Deep inelastic scattering 

The deep inelastic scattering (DIS) cross section is 
calculated using the GRV98 parton distribution func- 
tions [39| . Additionally, we have included the corrections 
in the small Q 2 region developed by Bodek and Yangjioj 
for both NEUT and NUANCE. The implementation of 
the model is slightly different in NEUT and NUANCE. 

In NEUT, the DIS contribution is included for 
hadronic invariant masses W > 1.3 GeV/c 2 . The pion 
multiplicity is additionally restricted to be greater than 
one for 1.3 < W < 2 GeV/c 2 to avoid double-counting 
sources of single pion production that are already in- 
cluded in the resonance portion of the simulation. The 
multi-hadron final states are simulated with two models 
in NEUT: a custom-made program [4l| for events with 
W between 1.3 and 2.0 GeV/c 2 and P YTHI A / JETSET 
[42[ for events with W larger than 2 GeV/c 2 . 
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A restriction on pion multiplicity is not enforced by 
NUANCE. Instead, the DIS contribution slowly increases 
for W values starting at 1.7 GeV and becomes the only 
source of neutrino interactions above W > 2 GeV. This is 
done to create a smooth transition between the resonance 
and DIS models and ensure continuity in distributions 
of kinematics and hadron multiplicity in the region of 
overlap. 



5. Intra-nuclear interactions 

Following production, the intra-nuclear interactions of 
mesons and nucleons are simulated using a cascade model 
in which the particles are traced until they escape from 
the nucleus. 

Although we only use kinematic information from the 
final state muon in this analysis, the simulation of intra- 
nuclear interactions is important since the pions/protons 
emitted from the nucleus can be mis-reconstructed as 
muons. 

The inelastic scattering, charge exchange and absorp- 
tion of pions in nuclei are simulated. For inelastic scat- 
tering and charge exchange interactions, the direction 
and momentum of pions are affected. In the scatter- 
ing amplitude, Pauli blocking is also taken into account. 
A more detailed description of the intra-nuclear interac- 
tion simulations in NUANCE and NEUT can be found 
elsewhere (9l Il5j. 

A 25% (30%) uncertainty in the overall pion absorp- 
tion (charge exchange) cross section is assumed based on 
comparisons to pion-carbon scattering data [43[ . The un- 
certainty in proton re-scattering is estimated to be 10%. 
Hence, we apply a 10% error on the number of proton 
tracks observed in SciBar. Additionally, we set a 20% un- 
certainty on the NC/CC ratio, estimated from the model 
dependence of the lepton-mass effect in the small Q 2 re- 
gion. 



C. The Expected Number of Neutrino Events 

Table U summarizes the parameter choices used in 
NEUT and NUANCE for the comparisons presented 
here. We chose these parameters since they are the de- 
fault parameter settings used in Super-K, K2K and T2K 
(NEUT) and MiniBooNE (NUANCE). The different pa- 
rameter values result in different neutrino cross section 
predictions between the two. 

With the SciBooNE neutrino beam exposure of 0.99 x 
10 20 protons on target, the expected number of events in 
the SciBooNE detector for each neutrino interaction is 
listed in Table |TTJ Because of the difference in parameter 
choices (Table IJ) , NEUT predicts a larger QE and single 
pion rate than NUANCE. The difference in QE rate is 
largely coming from the choice of k values, and the dif- 
ference in single pion rate can be largely accounted for 
by the difference in Ml" assumptions. The major source 



TABLE I. Parameters used for neutrino interaction simula- 
tion. 



Parameter 


NEUT 


NUANCE 


fF 


217 MeV 


220 MeV 


E B 


25 MeV 


34 MeV 


m1 e 


1.21 GeV 


1.23 GeV 


K 


1.00 


1.022 




1.21 GeV 


1.10 GeV 


M aoK 


1.0 GeV 


1.0 GeV 




(DIS) 


1.3 GeV 



of the factor of two larger multi-pion/DIS rate in NEUT 
compared to NUANCE is the difference of multi pion 
production simulation in the range 1.3 < W < 2.0 GeV; 
NEUT simulate these events as DIS, while NUANCE uses 
a resonant production model. These differences in cross 
section predictions between similar models with perfectly 
reasonable parameter choices further highlight the in- 
herent uncertainty in neutrino generator predictions and 
stress the importance of additional neutrino interaction 
measurements in this region. 

Figure [3] shows the expected number of total CC 
interactions as a function of neutrino energy. One can see 
that the cross section prediction from NEUT is about 10 
- 20% larger than that from NUANCE across the range 
of SciBooNE energies. 



^0.14 



x10 
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— NUANCE 




4.5 5 

E v (GeV) 
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0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



0.5 1 1.5 



2.5 



3.5 



4.5 5 
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FIG. 3. (color online). The neutrino energy distributions of 
CC interactions at the SciBar detector generated by NEUT 
(red dashed line) and NUANCE (black solid line). Top plot 
is the number of interactions per unit POT, and the bottom 
plot shows the ratio between the NEUT and NUANCE pre- 
dictions. 

The expected number of v„ CC interactions in the 
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TABLE II. The expected number and fraction of events in each neutrino interaction estimated by NEUT and NUANCE at the 
SciBooNE detector location with the neutrino beam exposure of 0.99 x 10 20 protons on target. The 10.6 ton fiducial volume 
of the SciBar detector is assumed. CC and NC interactions are abbreviated as CC and NC, respectively. 







NEUT 




NUANCE 


Interaction Type 


# Events 


Fraction(%) 


# Events 


Fraction(%) 


CC quasi-elastic 


53,038 


41.5 


47,573 


43.6 


CC single n via resonances 


29,452 


23.0 


25,863 


23.7 


CC coherent n 


1,760 


1.4 


1,736 


1.6 


CC multi-pion, DIS, etc 


6,834 


5.3 


3,140 


2.9 


NC total 


36,836 


28.8 


30,734 


28.2 


Total 


127,920 


100.0 


109,046 


100.0 



10.6 ton SciBar fiducial volume (Sec. HVB]) are 9.11 x 10 4 
and 7.83 x 10 4 for the NEUT and NUANCE predictions, 
respectively. 

IV. CC EVENT RECONSTRUCTION AND 
SELECTION 

To measure the rate of CC inclusive interactions, we 
use neutrino events occurring in SciBar with a muon in 
the final state. We select muons originating in the SciBar 
fiducial volume (FV), defined to be ±130 cm in both the 
x and y dimensions, and 2.62 < z < 157.2 cm, a total 
mass of 10.6 tons. 

We describe the reconstruction of muon tracks in 
Sec. IIV A[ the event selections and details the of sub- 
samples in Sec. IIV Bl and comparison of the data to the 
MC predictions in Sec. HVCl 



A. Track Reconstruction 

1. Track finding 

The first step of the event reconstruction is to search 
for two-dimensional tracks in each view of SciBar using 
a cellular automaton algorithm fjij . Three dimensional 
tracks are reconstructed by matching the timing and z- 
edges of the two dimensional projections; differences be- 
tween two two-dimensional projections are required to 
be less than 50 ns, and the z-edge difference must be less 
than 6.6 cm for both upstream and downstream edges. 
The two dimensional tracks in MRD arc independently 
reconstructed using hits in the MRD clustered within a 
50 ns timing window. Three dimensional tracks in the 
MRD are reconstructed by matching the timing of the 
two dimensional projections. 

Then, if the downstream edge of a SciBar track lies in 
last two layers of SciBar, we search for a track or hits in 
the MRD that are matched with the SciBar track. For 
matching an MRD track to a SciBar track, the upstream 



edge of the MRD track is required to be on either one 
of the first two layers of the MRD, and to be within 
30 cm of the projected entry point of the SciBar track 
into the MRD in each view. A SciBar track matched with 
MRD is defined as a SciBar-MRD matched track. The 
matching criteria impose a muon momentum threshold 
of 350 MeV/c. A more detailed description of the track 
reconstruction can be found elsewhere [15j. 



2. Track Classification 

We define three types of tracks used in this analy- 
sis: SciBar-stopped, MRD-stopped and MRD-penetrated 
tracks, as shown in Fig. 2J 



SciBar 


EC 


MRD 


i SciBar-FV 








1 SciBar-stopped i 






MRD-stopped 
MRD-penetrated 











FIG. 4. (color online). Types of tracks used for this CC 
interaction measurement. 

Tracks with both edges contained in the SciBar FV are 
classified as SciBar-stopped tracks. MRD-stopped and 
MRD-penetrated tracks are both subsets of the SciBar- 
MRD matched sample. An MRD-stopped track is se- 
lected by requiring the downstream edge of the track to 
be within \x\ < 132 cm , \y\ < 111 cm , and z < (last 
layer of the MRD). An MRD-penetrated track is selected 
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by requiring additional hits at the most downstream scin- 
tillator of the MRD. Tracks which exit from the side of 
the MRD are not used in this analysis. 

3. Particle Identification 

The SciBar detector has the capability to distinguish 
protons from other particles using dE/dx. We define a 
muon confidence level (MuCL) using the observed energy 
deposit per layer for all reconstructed tracks |15| . Tracks 
with MuCL greater than 0.05 are considered muon-like 
(or pion-like) and the others are classified as proton-likc. 

4- Muon Kinematics Reconstruction 

The slopes of the muon angles with respect to the beam 
in the two SciBar views are used to calculate the three 
dimensional muon angle with respect to the beam (0^). 

According to the MC simulation, about 30 % of SciBar- 
stopped tracks are expected to be backward-going. To 
identify these backward-going tracks, we use the de- 
layed timing signal produced by the decay electrons from 
stopped muons. A track is identified to be backward- 
going if it has (1) at least one delayed hit at t > 200 nsec 
at the upstream ends of the two dimensional views, and 
(2) no delayed hit at t > 200 nsec at the downstream 
ends. We impose this requirement for both top- and 
side- PMT signals from the track endpoints to remove 
hits due to random noise and after pulsing from the 
PMTs. The performance of this identification is esti- 
mated using a MC sample of CCQE events with 1 recon- 
structed track. In this sample, the efficiency of selecting 
backward-going tracks is ~ 57%, and the probability of 
mis-reconstructing a forward track as backward is ~1%. 
The loss of efficiency is predominantly due to decay elec- 
trons emitted at a large angle, producing hits in cither 
the top or side PMTs, but not both. 

This track direction identification is only applied to 
SciBar-stopped tracks. All tracks reaching the MRD 
are assumed to be forward-going, because the fraction 
of backward tracks is estimated to be small (~ 4%) and 
also a similar tagging of decay electrons is not possible 
in the MRD since the electrons stop and are undetected 
within the steel plates. 

The kinetic energy of the muon (Eki n ) is calculated by 
the range and the expected energy deposition per unit 
length (dE/dx) in the detector materials 

Ejg- — E SciBar 4- E EC + E Wal1 + E^ IRD (1) 

where E SciBar , E EC , E Walt and E MRD are the expected 
energy deposit by muons in SciBar, the EC, the wall of 
the dark box between the EC and MRD, and the MRD, 
respectively. 

For SciBar-stopped tracks, E EC , E Wal1 and E MRD 
are set to 0, and E SclBar is calculated by a range to 
energy look-up table based on the MC simulation. 



For MRD-stopped tracks, energy deposits in SciBar, 
EC and the wall are computed as E SctBar — 
2.04 MeV/cm x L SB (cm), E EC = 90.8/cos6» M MeV 
and E WaU = 3.3/ cos 9^ MeV, where L SB is the recon- 
structed track length in SciBar. E MRD is calculated by 
a range to energy look-up table based on the MC simu- 
lation. 

The average muon angular resolution is 0.9° for all 
samples. The muon momentum resolutions are 15 MeV/c 
for SciBar-stopped and 50 MeV/c for MRD-stopped 
tracks, respectively. 



B. Event Selection and Classification 

1. Event Selection 

We select the highest momentum track with 
MuCL > 0.05 in an event as a muon candidate. We 
also require the reconstructed momentum to be greater 
than 0.25 GeV/c to reject short proton or pion tracks 
from neutral current interactions. Then, we require the 
upstream edge of the muon candidate to be in the SciBar 
FV. 



2. Event Classification 

Events with muon track candidates are subdivided into 
three sub-samples: SciBar-stopped, MRD-stopped and 
MRD-penetrated samples, according to the track classi- 
fication given in Sec. II V A~2l for the muon candidate. 

a. SciBar-stopped sample The SciBar-stopped sam- 
ple provides the lowest energy sample; the mean en- 
ergy of neutrinos in this sample is 1.0 GeV. Accord- 
ing to the simulation, the purity of CC interactions 
in this sample is 85%. Impurities are due to NC 
interactions (~ 7%), interactions occurring in the sur- 
rounding material (~ 5%) and 77 M CC interactions (~ 
1.5%). Figure [5] shows the distributions of the recon- 
structed muon momentum (p^) and angle (# M ) for SciBar- 
stopped muons. The backward-going muon identifica- 
tion described in Sec. \TV A~4l is used for this sample. All 
tracks which are not identified as backward-going are ac- 
counted as forward-going tracks for both data and MC 
samples. The expected number of events in each in- 
teraction mode is summarized in Table Mil The con- 
tamination of the cosmic-ray backgrounds is estimated, 
using off-beam data, to be 671.2 ± 11.6 events, where 
the uncertainty comes from the statistics of the off-beam 
data sample. This rather large contamination (~ 5%) 
is due to cosmic-rays penetrate SciBar vertically, which 
leave sparse hits and are often mis-reconstructed as short 
tracks. This effect is negligible for the other sub-samples. 

b. MRD-stopped sample The MRD-stopped sample 
has the largest statistics among the three. The mean en- 
ergy of neutrinos in the MRD-stopped sample is 1.2 GeV. 
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FIG. 5. (color online). Distributions of reconstructed mo- 
mentum (top) and angle (bottom) of the muon candidate in 
the SciBar-stopped sample. The MC prediction is based on 
NEUT and absolutely normalized by the number of POT. 
The total and flux systematic errors on the MC predictions 
are separately shown. The dominant source of the total error 
is the cross section uncertainty. 



TABLE III. The expected number and fraction of events in 
each neutrino interaction type for the SciBar-stopped sample, 
as estimated by NEUT and NUANCE. The external events 
comprise interactions in the EC, MRD and surrounding ma- 
terial. 



Interaction 




NEUT 


NUANCE 


type 


Events 


Fraction(%) 


Events 


Fraction(%) 


CC QE 


5818 


47.4 


5511 


50.8 


CC res. lvr 


3694 


30.1 


3226 


29.8 


CC coh. Itt 


119 


1.0 


123 


1.1 


CC other 


914 


7.4 


350 


3.2 


NC 


916 


7.4 


842 


7.8 


All non-^ M 


188 


1.5 


161 


1.5 


External 


629 


5.1 


629 


5.8 


Total 


12278 




10,842 





According to the simulation, the purity of CC inter- 
action in this sample is 91%. Impurities are from neu- 
trino interactions in the EC/MRD which back-scatter 
(~ 4%), NC interactions (~ 3%) and 77 p CC inter- 
actions (~ 2%). Figure [6] shows the distributions of the 
reconstructed muon momentum (p M ) and angle (8^) for 
MRD-stopped muons. The expected number of events 
in each interaction mode is summarized in Table IIVI 
The contamination of the cosmic-ray backgrounds is es- 
timated to be 54.6 ± 3.3 events. 



|i 1800 
S1600 



DATA 

MC total error 
MC flux error 
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DATA 

MC total error 
I MC flux error 



80 100 120 140 160 180 
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FIG. 6. (color online). Distributions of reconstructed momen- 
tum (top) and angle (bottom) of the muon candidate in the 
MRD-stopped sample. The MC prediction is based on NEUT 
and absolutely normalized by the number of POT. The total 
and flux systematic errors on the MC predictions are sepa- 
rately shown. The dominant source of the total error is the 
cross section uncertainty. 

c. MRD -'penetrated sample The MRD-penetrated 
sample is the highest energy sample among the three. 
The mean energy of neutrinos is 2.4 GeV. Although the 
track angles can be measured, we do not have the capa- 
bility to reconstruct the muon momentum for the tracks 
which exit the MRD. However, this sample can provide 
the normalization for the highest energy region. Hence, 
this sample is also used for the neutrino interaction rate 
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TABLE IV. The expected number and fraction of events in 
each neutrino interaction type for the MRD-stopped sample, 
as estimated by NEUT and NUANCE. The external events 
comprise interactions in the EC, MRD and surrounding ma- 
terial. 



Interaction 




NEUT 


NUANCE 


type 


Events 


Fraction(%) 


Events 


Fraction(%) 


CC QE 


10341 


56.1 


8385 


52.3 


CC res. lvr 


4789 


26.0 


4839 


30.2 


CC coh. Itt 


659 


3.6 


633 


3.9 


CC other 


1010 


5.5 


600 


3.7 


NC 


577 


3.1 


569 


3.5 


All non-i/ M 


320 


1.7 


281 


1.8 


External 


729 


4.0 


729 


4.5 


Total 


18,427 




16,036 





TABLE V. The expected number and fraction of events in 
each neutrino interaction type for the MRD-penetrated sam- 
ple, as estimated by NEUT and NUANCE. The external 
events comprise interactions in the EC, MRD and surround- 
ing material. 



Interaction 




NEUT 


NUANCE 


type 


Events 


Fraction(%) 


Events 


Fraction(%) 


CC QE 


2428 


60.0 


1943 


57.0 


CC res. Itt 


1008 


24.9 


976 


28.6 


CC coh. Itt 


140 


3.5 


130 


3.8 


CC other 


356 


8.8 


255 


7.4 


NC 


1.5 


0.04 


2.3 


0.07 


All non-^ M 


89 


2.2 


75 


2.2 


External 


27 


0. 7 


27 


0.8 


Total 


4049 




3407 





measurement. According to the simulation, the purity 
of CC interaction in this sample is 97%. Impurities 
mostly come from 77 M CC interactions (~ 2%). Figure [7] 
shows the distributions of the reconstructed muon an- 
gle of the MRD-penetrated muons. The expected 
number of events in each interaction mode is summa- 
rized in Table |V] The contamination of the cosmic-ray 
backgrounds is estimated to be 16.6 ± 1.8 events. 
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from the NEUT based MC simulation. By combining 
these three samples, we can obtain fairly uniform accep- 
tance for neutrinos above 0.4 GcV. 
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FIG. 7. (color online). Distribution of reconstructed angle 
of the muon candidate in the MRD-penetrated sample. The 
MC prediction is based on NEUT and absolutely normalized 
by the number of POT. The total and flux systematic errors 
on the MC predictions are separately shown. 



FIG. 8. (color online). (Top) Number of CC events in the 
SciBar FV as a function of true E v , predicted by the NEUT 
based simulation. The number of selected events in each sub- 
sample are also shown. (Bottom) Detection efficiency as a 
function of true neutrino energy for each sub-sample. 



Efficiency Summary 



C. Data Comparison to the MC prediction 



Figure [5] shows the efficiency of CC events as a function 
of true neutrino energy for each sub-sample, estimated 



Table IVII shows the number of events obtained from 
data and the predictions from NEUT and NUANCE 
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based MC simulations. The contamination of cosmic- 
ray backgrounds is estimated using the off-beam data, 
and have been subtracted from the data. For the total 
number of events from the three sub-samples, we find 
a data/MC normalization factor of 1.08 for the NEUT 
prediction, and 1.23 for the NUANCE prediction. 

To compare the MC predictions with data, the neu- 
trino energy(Ev) and the square of the four-momentum 
transfer(Q 2 ) are the key variables since a flux variation 
is purely a function of E v while a variation of the cross 
section model typically changes the Q 2 distribution. We 
reconstruct these variables assuming CC-QE interaction 
kinematics. The reconstructed E v is calculated as 

Erec _ ™ 2 P ~ K ^ E B f - ml + 2(m» - E B )E ii 

where m p , m n and m M are the mass of proton, neutron 
and muon, respectively, E^ is the muon total energy, and 
E B is the nuclear potential energy. The reconstructed Q 2 
is given by, 

Q 2 rec = 2E: ec (E^ - Pll cos e„) - m\. (3) 

Figure |H] shows the distributions of E r v ec and Q 2 ec for 
the SciBar-stopped and MRD-stopped samples. In these 
plots, data points are compared with the NEUT and NU- 
ANCE based MC predictions. We find that the data are 
consistent with the MC predictions within the systematic 
uncertainties. 



V. CC INTERACTION RATE ANALYSIS 
A. Method 

To calculate the CC inclusive interaction rate and cross 
section versus energy, we re-weight the predictions of 
NEUT or NUANCE based simulations in true energy 
bins by factors that are found to give the best agreement 
with the kinematic distributions for data versus MC pre- 
diction. 

The p^ vs. (p u -9 n ) distributions from the SciBar- 
stopped and the MRD-stopped samples, and M distri- 
bution from the MRD-penetrated sample are simulta- 
neously used for this measurement. Figure [TU] shows 
the Pp-Ofj, distributions of the SciBar-stopped and MRD- 
stopped samples, while the 6* M distribution for the MRD- 
penetrated sample is shown in Fig.[7J Events in the same 
Pfj,-9fi bins but in different sub-samples are not summed 
together, but treated as separate p a -9^ bins in the anal- 
ysis, and only bins with at least 5 entries are used for the 
fit. The total number of p^-9^ bins is 159; 71 from the 
SciBar-stopped, 82 from the MRD-stopped and 6 from 
the MRD-penetrated samples. 

We define 6 rate normalization factors (fa,--- , /s) 
which represent the CC interaction rate normalized to 
the MC prediction for each true energy region defined in 
Table |VTTJ The events at E v < 0.25 GeV are not used 



since these events are below our detection efficiency as 
shown in Fig. [SJ and also the fraction of these low energy 
interactions are negligibly small (< 1%) at the BNB flux. 
We calculate these rate normalization factors by compar- 
ing the MC predictions to the measured CC interaction 
rate. For each energy region, we generate the MC tem- 
plates for the p u -9 u distributions in each event sample; 
n Pjed j g ^ e p rec jj c ^ e( j num ber of events in the j-th p^-9^ 
bin, corresponding to energy bin i. The expected number 
of events in each p^-9^ bin, Nj re , is calculated as 

E u bins 

N pre d= /,<"'. (4) 

i 

Figures QT] and [12] are MC templates of the p^-9^ 
distributions for the SciBar-stopped and MRD-stopped 
samples. We see that there is a large contribution 
in the SciBar-stopped sample of events with E v below 
0.75 GeV. Hence, this sample is essential to determine the 
rate normalization factors in the low energy regions. The 
Pfj,-8fi distributions of the MRD-stopped sample clearly 
depends on E V1 up to 1.75 GeV. However, most of the 
events in the MRD-stopped sample with E v > 1.75 GeV 
have small reconstructed p M . These are events with en- 
ergetic pion or proton tracks that are mis-reconstructed 
as muons. Due to the weak constraint from the MRD- 
stopped sample on events with E v > 1.75 GeV, the 
MRD-penetrated sample is included in the fit since about 
2/3 of the events in this sample have E u > 1.75 GeV as 
shown in Fig. [131 

We find the rate normalization factors (fo,-- - ,/s) 
which minimize the \ 2 value defined as: 

Nbins 

X 2 =Y, (Nfs-Nr d )(V sys + V stat )- k \N^-Nr d )- 

(5) 

Here, N?J%> and are the observed and predicted 

numbers of events in the j(k)-th p fl -9 u bin, and Njffc) 
is a function of the rate normalization factors as shown 
in Eq. V sys is the covariance matrix for systematic 
uncertainties in each Pu-9^ bin, and Vstat represents the 
statistical error. We have a total of 159 bins, so V sys and 
V s tat are 159 x 159 dimensional matrices. The details of 
evaluating V sys are described in the following section. 



B. Systematic Errors 

The sources of systematic error are divided into four 
categories: neutrino beam (i), neutrino interaction mod- 
els (ii), intra- nuclear interaction model (iii), and detec- 
tor response and neutrino interaction models outside of 
SciBar (iv) . Table I Villi shows the list of the systematic 
uncertainties considered in this analysis, which were de- 
scribed in the previous sections. 
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TABLE VI. The number of events in each sub-sample from the data and the predictions from NEUT/NUANCE-based MC. 
The numbers in parentheses show the ratio between the data and the predictions. The cosmic-ray backgrounds are estimated 
from off-timing data and subtracted from the data. 



Sample 


SciBar-stopped 


MRD-stopped 


MRD-penetrated 


Total 


Data 


13588.8 


20236.4 


3544.4 


37369.6 


NEUT 


12278.3(1.11) 


18426.3(1.10) 


4049.0(0.88) 


34753.6(1.08) 


NUANCE 


10841.9(1.25) 


16036.2(1.26) 


3407.5(1.04) 


30285.6(1.23) 



3000- 
2500^ 
2000^ 
1500- 
1000^ 
500^ 

oj^ 



DATA 
NEUT 
NUANCE 




0.5 1 



IS"—* 



ET(GeV) 





• DATA 
NEUT 
NUANCE 




EHGeV) 

FIG. 9. (color online). Top: E r v ec (left) and Ql ec (right) of the SciBar-stopped sample. Bottom: E r v ec (left) and Ql ec (right) 
of the MRD-stopped sample. The NEUT and NUANCE predictions are absolutely normalized by the number of POT. The 
filled regions show the systematic uncertainties on the MC predictions based on NEUT. The systematic uncertainty for the 
NUANCE prediction is similar to that of the NEUT prediction and not shown. 



The aim of this fit is to constrain the variations due to 
flux (i) and cross section (ii) by measuring the rate nor- 
malization factors as a function of E u , /j. Uncertainties 
from (i) and (ii) are factorized into two parts: variations 
which change the rate normalization of CC events in each 
E v region ( "normalization" ) and variations which change 
the Pp,-9^ distributions but not the normalization of CC 
events ("p^-Op shape"). The former, normalization, is 
corrected by the use of /j, so V sys corresponds only to 
Pn'O^i shape uncertainties for flux and cross section. 

The p^-O^ shape uncertainties are estimated by re- 



normalizing the variation for each E u region as follows. 
First, we generate a new prediction corresponding to a 
single systematic variation, n'^. Here, the prime denotes 
a systematic variation, with i and j representing as before 
the energy bin i and p^-O^ bin j. 

Then, the predicted event rate, iVj, corresponding to 
this systematic variation is: 

E u bins 

N j= E fi<i^ ( 6 ) 

i 

where Ri is the renormalization factor, which is the ratio 
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FIG. 10. Data distributions of p M vs. # M for the SciBar-stopped (left) and MRD-stopped (right) samples. The size of boxes is 
proportional to the number of entries. 



TABLE VII. Energy regions for the CC interaction rate measurement. These energy regions are in terms of the true neutrino 
energy from the MC. 



Parameter 


fo 


h 


h 


h 


h 


h 


E v range (GeV) 


0.25 - 0.5 


0.5 - 0.75 


0.75 - 1.0 


1.0 - 1.25 


1.25 - 1.75 


1.75+ 



of the total number of predicted events in each E v bin 
between central value and systematically varied predic- 
tions, defined as: 

n _ £±j _L (7\ 

The predictions for central value Nj and n,-j are the 

same as the Nf red and nfj ed used in Eq. (HI). 

j i j 1 * * 

Unlike the flux and cross section uncertainties, intra- 
nuclear interaction (iii) and detector uncertainties (iv) 
are independent of the CC interaction rate. For these 
sources, the uncertainties are simply calculated as: 

E w bins 

N'i= E f<- ( 8 ) 

i 

To generate the V sys matrix, the systematically varied 
rate, N-j, is compared to the central value event rate. 
To better estimate some uncertainties, many variation 
predictions ( "draws" ) are used, each corresponding to a 
unique set of underlying parameters. To estimate the er- 
ror from the Ma uncertainty, for example, we generate 
1000 different sets of MC expectations of riy, each cor- 
responding to the different value of Ma randomly drawn 
from the estimated uncertainty of Ma- Then, the error 



matrix representing the uncertainty from Ma is calcu- 
lated as: 

1 M 

^ = mD^~«"^ ( 9 ) 
i 

where k denotes the index of random draws and M de- 
notes the total number of draws. We estimate the er- 
rors from each source (i) - (iv) individually, then add the 
matrices together to generate a total covariance matrix, 
V sys . We assume the same size fractional error for both 
NEUT and NUANCE based predictions, and the total 
covariance matrices for them are calculated by scaling 
the fractional errors with the predictions. 

VI. RESULTS AND DISCUSSION 
A. Energy Dependent CC Interaction Rate 

1. Measurement of CC Interaction Rate 

Tables IIXI and [X] show the best fit parameters and 
their errors from the spectrum fit with the six energy de- 
pendent rate normalization factors for NEUT and NU- 
ANCE predictions, respectively. The ^/DOF, where 
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FIG. 11. The MC templates of p M vs. 6> M for the SciBar-stopped and MRD-stopped sample for the three lowest E v regions. 
The normalization factors are common between the SciBar-stopped and MRD-stopped samples. 
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FIG. 12. The MC templates of p M vs. 9^ for the SciBar-stopped and MRD-stopped sample for the three highest E v regions. 
The normalization factors are common between the SciBar-stopped and MRD-stopped samples. 
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TABLE VIII. List of systematic uncertainties considered. 



Category 


Error Source 


Variation 


Description 




7r + /ir~ production from p-Be interaction 


Spline fit to HARP data [19] 


Sec. |II B| 




K + / K° production from p-Be interaction 


Tables VIII and IX in Ref. [21] 


Sec. |II B| 


(i) 


Nucleon and pion interaction in Be/Al 


Table XIII in Ref. [21] 


Sec. |II B| 


Flux 


Horn current 


±1 kA 


Sec. |II B| 




Horn skin effect 


Horn skin depth, ±1.4 mm 


Sec. |II B| 




Number of POT 


±2% 


Sec.lHBj 




Fermi surface momentum of carbon nucleus 


±30 MeV 


Sec. HUB 11 




Binding energy of carbon nucleus 


±9 MeV 


Sec. HUB 11 


(") 


CC-QE M A 


±0.22 GeV 


Sec. HUB 11 


Neutrino 


CC-QE k 


±0.022 


Sec. HUB 11 


interaction 


CC-Itt M a 


±0.28 GeV 


Sec. HUB 21 




CC-Itt Q 2 shape 


Estimated from SciBooNE data 


Sec. HUB 21 




CC-coherent-7r Ma 


±0.28 GeV 


Sec. HUB 31 




CC-multi-7r M A 


±0.52 GeV 


Sec. HUB 41 




A re-interaction in nucleus 


±100 % 


Sec. HUB 21 


(iii) 


Pion charge exchange in nucleus 


±20 % 


Sec. HUB 51 


Intra-nuclear 


Pion absorption in nucleus 


±35 % 


Sec. HUB 51 


interaction 


Proton re-scattering in nucleus 


±10 % 


Sec. HUB 51 




NC/CC ratio 


±20 % 


Sec. HUB 51 




PMT 1 p.e. resolution 


±0.20 


Sec. |II D| 




Birk's constant 


±0.0023 cm/MeV 


Sec. |II D| 




PMT cross-talk 


±0.004 


Sec. |II D| 


(iv) 


Pion interaction cross section in the detector material 


±10 % 


Sec. |II D| 


Detector 


dE/dx uncertainty 


±3%(SciBar,MRD), ±10%(EC) 


Sec. |II D| 


response 


Density of SciBar 


±1 % 


Sec. |II C| 




Normalization of interaction rate at the EC/MRD 


±20 % 


Sec.|IIlA| 




Normalization of interaction rate at the surrounding materials 


±20 % 


Sec.|lIlA| 




Contamination of cosmic-ray backgrounds 


Estimated from off-beam data 


Sec.HVB2l 



DOF is the number of degree of freedom for NEUT 
and NUANCE predictions are respectively 225.1/159 and 
560.7/159 before fitting, and they are 161.2/153 and 
173.6/153 after fitting. We obtain reasonable \ 2 values 
after the fit. We also tried to tune Q 2 distributions by 
fitting cross section parameters such as Ma and /t, how- 
ever, these give similar \ 2 values. So we decided to make 
fits with simple factors for true E u regions, as described 
above, because they are less model dependent. 

Figure [T4l shows the distributions of 6^ E r v ec and 
Qrec of the MRD-stopped sample, after applying the rate 
normalization factors obtained in this analysis. We esti- 
mate the constrained systematic error for each distribu- 
tion in the same way as described in the Sec. lVBl We also 
propagate the errors of the scale factors (fi) to the dis- 
tributions. The errors on fi obtained from the fit include 
the shape error from all the flux and the cross section 
uncertainties, and the absolute error from all the intra- 
nuclear interaction and detector response uncertainties, 
as they are included into the error matrix (V sys ). The 
errors shown in these plots are the quadrature sum of 



TABLE IX. Spectrum fit result for the NEUT prediction. The 
first two rows are the best fit values and associated uncertain- 
ties. The subsequent columns and rows represent the corre- 
lation coefficients for each parameters. 





fo 


h 


h 


h 


u 


h 


Best fit 


1.042 


1.032 


1.234 


1.290 


1.193 


0.789 


Error 


0.186 


0.095 


0.064 


0.084 


0.104 


0.083 


fo 


1.0000 


0.2951 


0.1653 


-0.2224 


-0.3646 


0.0438 


h 


0.2951 


1.0000 


0.0591 


-0.3674 


-0.5311 


-0.2187 


h 


0.1653 


0.0591 


1.0000 


-0.2317 


-0.1511 


-0.3301 


fz 


-0.2224 


-0.3674 


-0.2317 


1.0000 


-0.0011 


0.2534 


u 


-0.3646 


-0.5311 


-0.1511 


-0.0011 


1.0000 


-0.0680 


h 


0.0438 


-0.2187 


-0.3301 


0.2534 


-0.0680 


1.0000 



those constrained systematic errors and errors of the fi. 
We find that both NEUT and NUANCE predictions well 
reproduce the data distributions within the errors of this 
analysis. Also, we confirm that the constraint by this 
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FIG. 13. (color online). The MC templates of M for the 
MRD-penetrated sample. The black line shows the total MC 
prediction. The predictions for E v > 1.75 GeV, 1.25 < E v < 
1.75 GeV and 1.0 < E v < 1.25 GeV are also shown. 



TABLE X. Spectrum fit result for the NUANCE prediction. 
The first two rows are the best fit values and associated un- 
certainties. The subsequent columns and rows represent the 
correlation coefficients for each parameters. 





fo 


fx 


h 


h 


/* 


h 


Best fit 


1.650 


1.312 


1.359 


1.378 


1.360 


0.898 


Error 


0.203 


0.093 


0.057 


0.075 


0.107 


0.085 


fo 


1.0000 


0.2260 


0.0716 


-0.2470 


-0.2916 


0.0019 


fi 


0.2260 


1.0000 


0.0244 


-0.4319 


-0.5275 


-0.2410 


h 


0.0716 


0.0244 


1.0000 


-0.2225 


-0.1150 


-0.3403 


h 


-0.2470 


-0.4319 


-0.2225 


1.0000 


0.0910 


0.2003 


f* 


-0.2916 


-0.5275 


-0.1150 


0.0910 


1.0000 


-0.0569 


h 


0.0019 


-0.2410 


-0.3403 


0.2003 


-0.0569 


1.0000 



measurement can reduce the systematic uncertainty in 
most regions, compared to the original errors. 

The CC interaction rate in i-th true E v region, TZi, is 
calculated as: 



n,, = 



■pred 



■Pi 



(10) 



where J\ff red is the number of selected events predicted 
by the MC simulation, Pj is the CC inclusive purity, and 
ti is the CC inclusive efficiency We evaluate the er- 
rors on Pi/ti with all systematic uncertainties listed in 
Table I Villi Then we take the quadrature sum of the 
errors of fi and Pj/ej to estimate the total systematic 
error. Figure [TS] shows the obtained neutrino interac- 



tion rate normalized to the NEUT and NUANCE pre- 
dictions, with the full systematic error. The original flux 
and cross-section uncertainties are also shown in the plot. 
The numerical values of the interaction rate normalized 
to the MC predictions and its full errors are shown in 
Table IXII The uncertainty is about 6% for the rate at 
0.75 < E u < 1.0 GeV, where the CC interaction rate is 
maximum, and is about 15% for the lowest energy region. 

TABLE XI. CC inclusive interaction rate normalization 
factors to NEUT and NUANCE predictions. The size of the 
full systematic errors are also shown. 

Energy region CC rate normalization factor 

(GeV) NEUT NUANCE 



0.25 - 0.50 
0.50 - 0.75 
0.75 - 1.00 
1.00 - 1.25 
1.25 - 1.75 
1.75 - 



1.04 ±0.20 
1.03 ± 0.11 
1.23 ± 0.08 
1.29 ±0.10 
1.19 ± 0.11 
0.79 ± 0.08 



1.65 ±0.22 
1.31 ±0.11 
1.36 ±0.08 
1.38 ±0.09 
1.36 ±0.12 
0.90 ±0.09 



2. CC Inclusive Cross Section 

The major difference between the rate normalization 
factors NEUT and NUANCE is (a) the difference in the 
predicted neutrino interaction cross sections of the two 
models. This effect is compensated by the rate normal- 
ization factor obtained by this analysis, and thus has 
little effect on the absolute rates. 

However, there is also a second-order contribution to 
this difference from (b) the difference in the prediction of 
final state particle production and kinematics. This af- 
fects events with pion or proton tracks mis-reconstructed 
as muons, which ultimately affects the purity (Pi) and 
efficiency (e;) in Eq. (fl"0|). 

To estimate the effect of (b), we extract the absolute 
CC interaction cross section from the obtained rate nor- 
malization factors. Here, the difference between NEUT 
and NUANCE due to (a) is canceled, and the remaining 
differences are, in principle, due to the source (b). The 
cross section per nucleon on polystyrene target (CsHg) 
at each energy region is calculated as 



cc 



•pred 



P 



£i ■ T ■ $ , 



11) 



where i is the index of the energy regions used for the 
spectrum fit (see Table IVlT|) . < >i is the predicted 

flux averaged CC interaction cross section per nucleon, 
j^j-pred j g num ber of selected events predicted by the 
MC simulation, Pi is the CC inclusive purity, ti is the 
CC inclusive efficiency, Ti is the number of nucleons in 
the SciBar fiducial volume, and <J>i is the muon neutrino 
flux per unit area. 
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FIG. 14. (color online). Top: Reconstructed p M (left) and # M (right) of the MRD-stopped sample. Bottom: _E^ ec (left) and 
Qrec (right) °f the MRD-stopped sample. The NEUT and NUANCE predictions after the fit are shown, and the predictions 
are absolutely normalized by the number of POT. The filled regions show the systematic uncertainties of MC predictions based 
on NEUT. The area within dashed lines shows the NEUT predictions with their systematic uncertainties before the fit. The 
systematic uncertainty for the NUANCE prediction is similar to that of the NEUT prediction and not shown. 



Figure [16] show the extracted cross sections plotted 
with the original predictions from NEUT and NUANCE. 
In addition to the errors on Pj/e, as estimated for the 
rate measurements, we also estimate the errors on 
from the category (i) in the table. In the plot, we sep- 
arately show the errors of fi and the quadrature sum of 
fi, Pi, ej and $i errors. We confirm that the differences 
of the extracted CC interaction cross sections between 
NEUT and NUANCE are within the errors of fi. There- 
fore, the effect of source (b) is small and covered by the 
systematic uncertainty. The difference of the rate nor- 
malization factors is mostly caused by the cross section 
difference itself (source (a)). 



The obtained cross section values and their errors are 
summarized in Table IXIII The uncertainty on the cross 
section is about 10% at 0.75 < E v < 1.0 GeV, where the 
CC interaction rate is maximum, and is about 30% for 
the lowest energy region. 



B. Flux Integrated CC Interaction Rate 

1. Measurement of CC Interaction Rate 

We perform a fit with the method described in Sec. |V] 
using one rate normalization factor (ftot) which spans the 
entire energy region, instead of six regions. This is moti- 
vated by the SciBooNE-MiniBooNE joint disappear- 
ance analysis, in which we also evaluate the sensitivity 
using one rate normalization factor. 

The obtained rate normalization factors are 1.12 ± 
0.02 ± 0.03, and 1.29 ± 0.02 ± 0.03 for the NEUT and 
NUANCE predictions, respectively. The first errors show 
the errors from the fit, and the second show the errors 
from efficiency and purity uncertainties. Here, higher 
rate normalization factors compared to the simple ratio 
of the number of events (Table IVII) are obtained. This 
is because the "shape" systematic uncertainty fSec. lVBj) 
is calculated over the entire energy region, and it typi- 
cally becomes smallest where the CC interaction rate is 
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TABLE XII. Energy dependent CC inclusive cross section per nucleon on a polystyrene target (CsHs). Results based on 
NEUT and NUANCE based predictions are separately shown. 



Energy region 


Mean Energy 


Total Vfj, flux 


CC inclusive cross 


section (cm 2 /nucleon) 




(GeV) 


(GeV) 


(i/^/cm ) 


IN hi U 1 based 


IN U A1N Ljhj based 




0.25 - 0.50 


0.38 


(4.31 ±0.81) x 10 11 


(2.76 ± 0.75) x 10 -39 


(3.40 ± 0.96) x 10" 


39 


0.50 - 0.75 


0.62 


(5.09 ± 0.37) x 10 11 


(5.80 ± 0.75) x 10~ 39 


(6.39 ±0.81) x 1(T 


39 


0.75 - 1.00 


0.87 


(4.18 ±0.26) x 10 11 


(1.03 ±0.10) x 10 -38 


(1.01 ±0.09) x 10" 


38 


1.00 - 1.25 


1.11 


(2.63 ± 0.23) x 10 11 


(1.38 ±0.17) x 10~ 38 


(1.29 ±0.15) x 10" 


38 


1.25 - 1.75 


1.43 


(1.90 ± 0.27) x 10 11 


(1.62 ±0.29) x 10~ 38 


(1.56 ±0.28) x 10" 


38 


1.75 - 


2.47 


(0.62 ± 0.12) x 10 11 


(1.74 ±0.38) x 10 -38 


(1.66 ±0.37) x 10" 


38 



x10 



■39 




FIG. 15. (color online). CC interaction rate normalized to the 
NEUT and NUANCE predictions, obtained by the spectrum 
fit. The error bars show the full systematic uncertainties. The 
original flux and cross section errors are also shown separately. 



maximum. Hence, the obtained rate normalization fac- 
tors tend to be close to the results for 0.75 < E v < 
1.0 GeV in the energy dependent analysis (Sec. IVI Al) . 
The X 2 /DOF, for NEUT and NUANCE predictions are 
respectively 208.6/159 and 481.5/159 before fitting, and 
they are 173.0/158 and 183.0/158 after fitting. 



2. CC Inclusive Cross Section 

We extract the total CC inclusive cross section for the 
flux at E v > 0.25 GeV using a similar formula to Eq. (ITT1) : 



r , pred . 

C = Jtof < & CC > = 



ftofN Pred -P 



(12) 



The purity (P) and the efficiency (e) are estimated for 
the sum of SciBar-stopped, MRD-stopped and MRD- 
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FIG. 16. (color online). CC inclusive interaction cross section 
per nucleon on a polystyrene target (CsHs). The smaller error 
bars show the uncertainties of the rate normalization factors, 
and the larger error bars represents the total error including 
the flux uncertainties. 



penetrated samples, and obtained to be P = 90.0(89.1)% 
and e = 34.5(34.5)% with the NEUT(NUANCE) genera- 
tors. The integrated flux at E v > 0.25 GeV averaged 
over the SciBar FV for 9.9E19 POT is estimated to be 
<t> = (1.87 ± 0.14) x 10 12 (y M /cm 2 ), with mean energy of 
0.83 GeV. 

Using these purity, efficiency and flux predictions, 
we obtain the total CC cross sections per nucleon 
on a polystyrene target for a muon neutrino beam 
with mean energy of 0.83 GeV to be (8.51 ± 0.72) x 
10~ 39 cm 2 /nucleon and (8.45±0.71) x 10 -39 cm 2 /nucleon 
for NEUT and NUANCE based predictions, respectively. 
The error contains all systematic uncertainties from the 
purity and efficiency variation and absolute flux uncer- 
tainty. The error is dominated by the uncertainty on the 
total flux prediction (7.6%). We obtain consistent CC in- 
teraction cross sections from the two neutrino generator 
simulations. 



20 



TABLE XIII. Extracted number of total CC interaction and absolute cross sections per nucleon on a polystyrene target, based 
on the number of events in each sub-sample. The MC predictions are based on NEUT. 





Events 


Mean energy of 


Efficiency 


Purity 


Total CC events 


Flux integrated CC 


sample 


(Af obs ) 


selected sample [GeV] 


(c) 


(P) 


(Af obs P/e) 


cross-section (a [cm 2 /nucleon]) 


MRD-match 


28409.6 


1.43 


27.2% 


92.9% 


(9.71 ± 0.69) x 10 4 


(8.14 ±0.85) x 10" 39 


MRD-stop 


20236.4 


1.20 


18.5% 


91.2% 


(10.0 ±0.73) x 10 4 


(8.34 ±0.88) x 10~ 39 


MRD-pene 


3544.4 


2.42 


4.3% 


97.1% 


(7.94 ± 1.54) x 10 4 


(6.65 ± 1.29) x 10~ 39 



CC Inclusive Cross Section for Previous SciBooNE 
Results 



In previous SciBooNE publications, the total CC inter- 
action cross sections were used to normalize the results. 
We have used the number of MRD-matched and MRD- 
penetrated samples to extract the total CC cross sec- 
tion for CC coherent pion production measurements flij . 
and used MRD-stopped sample for NC neutral pion pro- 
duction measurements [H, [lj} ■ We extract the absolute 
CC inclusive cross section from the number of each sub- 
sample as: 



Af obs ■ P 



(13) 



where Af obs is the observed number of events in each 
sub-sample, P is the purity of CC interaction in the 
sample. 

Table IXlLTl is a summary of the absolute CC cross sec- 
tion extracted from the number of events in each sub- 
sample. Since we made minor improvements to the MC 
prediction and the reconstruction algorithms, the num- 
ber of events in each sub-sample are slightly different 
from the previous published results. The difference in 
the number of events are included in the systematic un- 
certainties. 

For the cross section values extracted from the MRD- 
matched and the MRD-stopped samples, the dominant 
source of uncertainty is the total flux ($) uncertainty 
(7.6%). The second largest uncertainty is the efficiency 
(e) uncertainty due to the cross section models. They 
are 4.0% and 5.0% in the MRD-matched and the MRD- 
stopped samples, respectively. 

For the MRD-penetrated sample, the dominant sys- 
tematic error from the efficiency variation is due to the 
flux model uncertainties. Approximately 40% of the 
MRD-penetrated sample is from neutrinos from kaon de- 
cay, which have a large uncertainty, as shown in Fig. [5] 
The error from the efficiency variation due to the flux 
model is estimated to be 18.1%. The second largest er- 
ror comes from the the total flux ($) uncertainty (7.6%). 

With these cross section values, one can convert the 
cross section ratios of our previous results to absolute 
cross sections. 



VII. CONCLUSIONS 

In conclusion, we have isolated three CC data sam- 
ples, measured and and extracted the CC inter- 
action rates and cross sections using flux and neutrino 
interaction simulations. 

We extract the CC interaction rates by fitting muon 
kinematics, with precision of 6-15% for the energy depen- 
dent and 3% for the energy integrated analyses. We con- 
firm that the distributions after fitting well reproduce the 
observed distributions with both NEUT and NUANCE 
based simulations. This result will be used to constrain 
the neutrino interaction rate for a SciBooNE-MiniBooNE 
joint disappearance analysis [13j . 

We also evaluate CC inclusive interaction cross sec- 
tions, and the results are consistent with both NEUT 
and NUANCE predictions. This confirms that the dif- 
ference in the observed rates normalized to the NEUT 
and NUANCE based predictions is mainly due to the 
cross section difference of the two simulators. The pre- 
cisions of the obtained cross sections are 10-30% for en- 
ergy dependent and 8% for the energy integrated anal- 
yses. This is the first measurement of the CC inclusive 
cross section on carbon around 1 GeV. The total CC in- 
teraction cross section values, as defined in our previous 
publications [l5T - [l7| . are also extracted. These cross sec- 
tion values can be used to convert previous measurements 
of cross section ratios of exclusive channels at SciBooNE 
to the absolute scale. The results may also be used to 
tune neutrino interaction models in the ~1 GeV region, 
which is relevant to various ongoing and future neutrino 
oscillation experiments. 



VIII. ACKNOWLEDGEMENTS 

We acknowledge the Physics Department at Chonnam 
National University, Dongshin University, and Seoul Na- 
tional University for the loan of parts used in SciBar and 
the help in the assembly of SciBar. We wish to thank 
the Physics Departments at the University of Rochester 
and Kansas State University for the loan of Hamamatsu 
PMTs used in the MRD. We gratefully acknowledge sup- 
port from Fermilab as well as various grants, contracts 
and fellowships from the MEXT and JSPS (Japan), the 
INFN (Italy), the Ministry of Science and Innovation 
and CSIC (Spain), the STFC (UK), and the DOE and 



21 



NSF (USA). This work was supported by MEXT and 
JSPS with the Grant-in-Aid for Scientific Research A 
19204026, Young Scientists S 20674004, Young Scien- 
tists B 18740145, Scientific Research on Priority Areas 
"New Developments of Flavor Physics" , and the global 



COE program "The Next Generation of Physics, Spun 
from Universality and Emergence" . The project was sup- 
ported by the Japan/U.S. Cooperation Program in the 
field of High Energy Physics and by JSPS and NSF under 
the Japan-U.S. Cooperative Science Program. 



J. Phys. Conf. Ser. 136, 042017 (2008) 




[9] D. Casper, 


Nucl. Phys. Proc. Suppl. 112, 161 (2002) 


arXiv:hep-ph/0208030 



Q. Wu et al. (NOMAD), Ph ys. Lett. B660, 19 (2008) 
arXiv:0711.1183 [hep-ex] | 

P. Adamson et al. (MINOS) 

Phys. Rev. D 81, 072002 (2010)| 
arXiv:0910.2201 [hep-ex] 



S. J. Barish et al, Phys. Rev. D 19, 2521 (1979) 
N. J. Baker et a/., |Phys. Rev. D 25, 617 (1982)| 

A. A. Agui lar-Arevalo et al. (Mini 

BooNE), |Phys. Rev. Lett 98, 231801 (2007) 



arXiv:0704.1500 [hep-ex 
Y. Itow et al. (T2K), | arXiv:hep-ex/0106019 
Y. Hayato, [Nucl. Phys 



Proc. 



G. 



Mitsuka 



Suppl. 112, 171 (2002)| 
(Super-Kamiokande) , 



al. 



(Kamiokande) 



Y Fukuda et 

|Phys. Lett. B335, 237 (1994)| 
Y. Ashie et al. (Super-Kamiokande), 
|Phys. Rev. D 71, 112005 (2005)||arXiv:hep-ex/0501064| 
M. H. Ahn et al. (K2K), 

|Phys. Rev. D 74, 072003 (2006)| |arXiv:hep-ex/06 06032 
Y. Nakajima (SciBooNE), |arXiv: 1010.5721 [hep-exj| 
A. A. Aguilar-Arevalo et al. (Mini- 

BooNE), |Phys. Rev. Lett. 103, 061802 (200971 

|arXiv:0903.2465 [hep-exj] ~ 
K. Hiraide et 

|Phys. Rev. D 78, 112004 (2008)| 
Y. Kurixnoto et 

Phys. Rev. D 81, 033004 (2010)| 



al. 



al. 



(SciBooNE), 
(SciBooNE), 



arXiv:0910.5768 [hep-ex] 
Y. Kurimoto 
Phys. Rev. D 81, 111102" 
arXiv:1005.0059 [hep-ex] 



et 

(2010)| 



al. 



(SciBooNE), 



S. Agostinclli 


et al. 


(GEANT4), 


Nucl. Instrum. Meth. A506, 250 (2003) 


M. G. Catanesi et al., 


Eur. Phys. J. C52, 29 (2007) 


arXiv:hep-ex/0702024| 


et al 


(E910), 


I. Chemakin 



Phys. Re v. C 77, 
arXiv:0707.2375 [nucl-exj] ~ 



A. A. 
BooNE), 



arXiv:0806.1449 [hep-ex] 



Aguilar-Arevalo et al. (Mini- 
Phys. Rev. D 79, 072002 (2009)] 



[22] 

[23] 
[24] 
[25] 



[26] 
[27] 

[28] 
[29] 



[30] 
[31] 

[32] 



[33] 
[34] 
[35] 
[36] 

[37] 



K. Nitta et aZ.,|Nucl. Instrum. Meth. A535, 147 (2004) 
|arXiv:hep-ex/0406023p 

H. Maesaka, Ph. D. thesis, Kyoto University (2005). 
M. Yoshida et aZ., |IEEE Trans. Nucl. Sci. 51, 3043 (2004) 
A. Heikkinen, N. Stepanov, and J. P. Wellisch, in Pro- 
ceedings of 2003 Conference for Computing in High- 
Energy and Nuclear Physics ( CHEP 03), La Jolla, 
California, 24-28 Ma r 2003 (2003) p. MOMT008, 
|arXiv:nucl-th/0306008| 

M. Hasegawa, Ph.D. thesis, Kyoto University (2006). 
R. A. Smith and E. J. Moniz, 

|Nucl. Phys. B43, 605 (1972)] 

E. J. Moniz et fll., |Phys. Rev. Lett. 26, 445 (1971) 
A. A. Aguilar-Arevalo et al. 
BooNE), |Phys. Rev. Lett. 100, 032301 (2008) 



(Mis 



arXiv:0706.0926 [hep-ex 
H 



and J. Arrington, 



S. Budd, A. Bodek, 
|arXiv:hep-ex/0308005 1 

X X Aguilar- Arevalo et al. (Mini- 

BooNE), |Phys. Rev. D 81, 092005 (2010)| 



i^rXiv : 1002.2680 [hep^ex] 

X Bodek, S~! A~vvakumov, R. Bradford, and 
H. S. Budd, |J. Phys Conf. Ser. 110, 082004 (2008)| 
|arXiv:0709.3538 [hep-exj] 



D. Rein and L. M. Sehgal,|Ann. Phys. 133, 79 ( 1981) 
D. Rein, |Z. Phys. C35, 43 (1987JP 



D. Rein and L. M. Sehgal JNucl. Phys. B 223, 29 (1983) 
D. Rein and L. M . Sehgal, Phys. Lett. B657, 207 (2007) 
|arXiv:hep-ph/0606185| 

K. Nakamura et al. (Particle Data Group) 
|J Phys. G37, 075021 (2010)] 

X X Aguilar-Arevalo et al. (MiniBooNE) 



[39] 

[40] 
[41] 

[42] 
[43] 



[44] 



|Phys. Lett. B664, 41 (2008)] |arXiv:0803.3423 [hep^ex 
M. Gluck, E. Reya, and A. Vo gt, 

|Eur. Phys. J. C5, 461 (1998)] |arXiv:hep -ph/9806 404] 
A. Bodek and U. K. Yang, |arXiv:hep-ex/0308007| 
M. Nakahata et al. (Kamiokande), 

|J. Phys. Soc. Jap. 55, 3786 (1986) | 

T. Sjostrand, |Comput. Phys. Commun. 82, 74 (1994)| 
D. Ashery et al, |Phys. Rev. C 23, 2173 (1981)| R. D. 
Ransome et al., |Phys. Rev. C 45, 509 (1992JP M. K. 
Jones et al., |Phys. Rev. C 48, 2800 (1993)| ~ 
J. L. Alcaraz-Aunion, Ph.D. thesis, Universidad Au- 
tonoma de Barcelona (2010). 



